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Abstract. We consider Glauber dynamics for the Ising model on the 
complete graph on n vertices, known as the Curie-Weiss model. It is 
well-known that the mixing-time in the high temperature regime (/3 < 1) 
has order n log n, whereas the mixing-time in the case f3 > 1 is exponen- 
tial in n. Recently, Levin, Luczak and Peres proved that for any fixed 
(3 < 1 there is cutoff at time 2 ^_^ n\ogn with a window of order n, 
whereas the mixing-time at the critical temperature /? = 1 is 0(n 3/ ' 2 ). 
It is natural to ask how the mixing-time transitions from G(nlogn) to 
0(n 3 / 2 ) and finally to exp(0(n)). That is, how does the mixing-time 
behave when (3 — /3(n) is allowed to tend to I as n — > co. 

In this work, we obtain a complete characterization of the mixing- 
time of the dynamics as a function of the temperature, as it approaches 
its critical point f3 c — 1. In particular, we find a scaling window of 
order 1/y/n around the critical temperature. In the high temperature 
regime, /3 = 1 — 5 for some < 8 < 1 so that S 2 n — > oo with n, the 
mixing-time has order (n/<5) log(<5 2 n), and exhibits cutoff with constant 
| and window size n/8. In the critical window, (3 = 1 ± 5 where S 2 n 
is O(l), there is no cutoff, and the mixing-time has order n 3 ^ 2 . At low 
temperature, /3 = l + <5for5>0 with 5 2 n — ► oo and S = o(l), there is 
no cutoff, and the mixing time has order j exp ((| + o(l))S 2 n). 



The Ising Model on a finite graph G = (V,E) with parameter (3 > 
and no external magnetic field is defined as follows. Its set of possible 
configurations is = {1,-1}^, where each configuration a G assigns 
positive or negatives spins to the vertices of the graph. The probability that 
the system is at a given configuration a is given by the Gibbs distribution 



where Z{j3) (the partition function) serves as a normalizing constant. The 
parameter (3 represents the inverse temperature: the higher (5 is (the lower 
the temperature is), the more [Xq favors configurations where neighboring 
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Cutoff only occurs at high temperature. 

spins are aligned. At the extreme case = (infinite temperature), the 
spins are totally independent and [iq is uniform over Q. 

The Curie- Weiss model corresponds to the case where the underlying 
geometry is the complete graph on n vertices. The study of this model (see, 
e.g., [7], [8], [9], [11]) is motivated by the fact that its behavior approximates 
that of the Ising model on high-dimensional tori. It is convenient in this 
case to re-scale the parameter (3, so that the stationary measure \i n satisfies 



The heat-bath Glauber dynamics for the distribution fi n is the following 
Markov Chain, denoted by (Xf). Its state space is O, and at each step, a 
vertex x S V is chosen uniformly at random, and its spin is updated as 
follows. The new spin of x is randomly chosen according to fj, n conditioned 
on the spins of all the other vertices. It can easily be shown that (Xf) is an 
aperiodic irreducible chain, which is reversible with respect to the stationary 
distribution \i n . 

We require several definitions in order to describe the mixing-time of the 
chain {X-t). For any two distributions <f>,ip on O, the total-variation distance 
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of (j) and if) is defined to be 

U - V||tv := sup \<f>(A) - ip(A)\ = IJ2 M<t) - • 

The (worst-case) total-variation distance of (Xt) to stationarity at time t is 
d n (t) := max \\P a (X t £ •) - ^ n || T v , 

where P CT denotes the probability given that Xq = a. The total-variation 
mixing-time of (Xt), denoted by t Mlx (e) for < e < 1, is defined to be 

W(e) := min{t : d n (t) < e} . 

A related notion is the spectral-gap of the chain, gap := 1 — A, where A is the 
largest absolute-value of all nontrivial eigenvalues of the transition kernel. 

Consider an infinite family of chains (X^ ) , each with its corresponding 
worst-distance from stationarity d n (t), its mixing-times t[^ x > etc. We say 
that (X^: ) exhibits cutoff iff for some sequence w n = o(t^x(|)) we have 
the following: for any < e < 1 there exists some c e > 0, such that 

4i"x( £ ) - - e) < c e w n for all n . (1.2) 

That is, there is a sharp transition in the convergence of the given chains to 
equilibrium at time (1 + o(l))t^x(i)- I n this case, the sequence w n is called 
a cutoff window, and the sequence 4«x(i) i s called a cutoff point. 

It is well known that for any fixed (3 > 1, the Glauber dynamics (Xt) 
mixes in exponential time (cf., e.g., [10]), whereas for any fixed j3 < 1 
(high temperature) the mixing time has order nlogn (see [2] and also [3]). 
Recently, Levin, Luczak and Peres [11] established that the mixing-time at 
the critical point (3 = 1 has order n 3 / 2 , and that for fixed < /3 < 1 there is 
cutoff at time 2 , 1 1 _^ nlog n with window n. It is therefore natural to ask how 
the phase transition between these states occurs around the critical (3 C = 1: 
abrupt mixing at time ( 2 (i_m + o(l))nlogn changes to a mixing-time of 
0(n 3 / 2 ) steps, and finally to exponentially slow mixing. 

In this work, we determine this phase transition, and characterize the 
mixing-time of the dynamics as a function of the parameter (3, as it ap- 
proaches its critical value C = 1 both from below and from above. The 
scaling window around the critical temperature (3 C has order l/^/n, as for- 
mulated by the following theorems, and illustrated in Figure 1. 

Theorem 1 (Subcritical regime). Let 5 = 5(n) > be such that 5 2 n — ► oo 
with n. The Glauber dynamics for the mean-field Ising model with parameter 
(3 = 1 — 5 exhibits cutoff at time \(n/ 5)\og(5 2 n) with window size n/S. In 
addition, the spectral gap of the dynamics in this regime is (1 + o(l))5/n, 
where the o(l)-term tends to as n — > oo. 



1 



JIAN DING, EYAL LUBETZKY AND YUVAL PERES 



Theorem 2 (Critical window). Let 5 = 5{n) satisfy 5 = 0(l/y/n). The 
mixing time of the Glauber dynamics for the mean-field Ising model with 
parameter (3 = 1±<5 has order n 3 / 2 , and does not exhibit cutoff. In addition, 
the spectral gap of the dynamics in this regime has order n~ 3 / 2 . 

Theorem 3 (Supercritical regime). Let 5 = 5{n) > be such that 5 2 n — ► oo 
with n. The mixing-time of the Glauber dynamics for the mean-field Ising 
model with parameter (3 = 1 + 5 does not exhibit cutoff, and has order 

n (n /< fl + g(x)\ , \ 
t cxp (n):=-exp^H ^{l^gjxj J **) ' 

where g[x) := (tanh(/3x) — x) / (1 — x tanh(/3x)) ; and £ is the unique positive 
root of g. In particular, in the special case 5^0, the order of the mixing 
time is j exp ((f + o{l))5 2 n), where the o(l)-term tends to as n — » oo. In 
addition, the spectral gap of the dynamics in this regime has order l/t exp (n). 

As we further explain in Section 2, the key element in the proofs of the 
above theorems is understanding the behavior of the sum of all spins (known 
as the magnetization chain) at different temperatures. This function of the 
dynamics turns out to be an ergodic Markov chain as well, and namely a 
birth- and- death chain (a one- dimensional chain, where only moves between 
neighboring positions are permitted). In fact, the reason for the exponential 
mixing at low-temperature is essentially that this magnetization chain has 
two centers of mass, +C,n (where C is as defined in Theorem 3), with an 
exponential commute time between them. Figure 2 demonstrates how the 
single center of mass around that this chain (rescaled) has at high near- 
critical temperature proceeds to split into two symmetric centers of mass 
that drift further and further apart as the temperature decreases. 

In light of this, a natural question that rises is whether the above men- 
tioned bottleneck between the two centers of mass at ±£n is the only reason 
for the exponential mixing-time at low temperatures. Indeed, as shown in 
[11] for the strictly supercritical regime, (3 > 1 fixed, if one restricts the 
Glauber dynamics to non-negative magnetization (known as the censored 
dynamics), the mixing-time becomes 0(nlogn) just like in the subcritical 
regime. Formally, the censored dynamics is defined as follows: at each step, 
a new state a is generated according to the original rule of the Glauber 
dynamics, and if a negative magnetization is reached (S(o~) < 0) then a is 
replaced by —a. Interestingly, this simple modification suffices to boost the 
mixing-time back to order nlogn, just as in the high temperature case, and 
thus raises the question of whether the symmetry between the high temper- 
ature regime and the low temperature censored regime applies also to the 
existence of cutoff. 
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Figure 2. The stationary distribution of the normalized 
magnetization chain (average of all spins) for the dynamics 
on n = 500 vertices. The center of mass at high tempera- 
tures (see j3 = 0.95) is at 0. Low temperatures feature two 
centers of mass at ±£ (where C, is the unique positive solution 
of tanh(/3x) = x), leading to the exponential mixing time. 

In a companion paper [6], we strengthen the result of [11] by showing that 
the scaling window of l/\/n exists also for the censored low temperature 
case, beyond which cutoff indeed occurs (yet at a different location than in 
the symmetric high temperature point). 

Theorem 4. Let 5 > be such that 8 2 n — > oo arbitrarily slowly with n. 
Then the censored Glauber dynamics for the mean field Ising model with 
parameter (3 = 1 + 8 has a cutoff at 

( 1 1 \ n . , r2 N 

tn= [2 + 2(ep/5-l)) 5 l ° g{6 H) 

with a window of order n/5. In the special case of the dynamics started from 
the all-plus configuration, the cutoff constant is [2(£ 2 /?/# — 1)] _1 (the order 
of the cutoff point and the window size remain the same). 

Theorem 5. Let 5 > be such that S 2 n — > oo arbitrarily slowly with n. 
Then the censored Glauber dynamics for the mean field Ising model with 
parameter 0=1 + 8 has a spectral gap of order 8/n. 

Recalling Theorem 1, the above confirms that there is a symmetric scaling 
window of order 1 / y/n around the critical temperature, beyond which there 
is cutoff both at high and at low temperatures, with the same order of 
mixing-time (yet with a different constant), cutoff window and spectral gap. 
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The rest of this paper is organized as follows. Section 2 contains a brief 
outline of the proofs of the main theorems. Several preliminary facts on the 
Curie- Weiss model and on one-dimensional chains appear in Section 3. Sec- 
tions 4, 5 and 6 address the high temperature regime (Theorem 1), critical 
temperature regime (Theorem 2) and low temperature regime (Theorem 3) 
respectively. 

2. Outline of proof 

In what follows, we present a sketch of the main ideas and arguments 
used in the proofs of the main theorems. We note that the analysis of the 
critical window relies on arguments similar to those used for the subcritical 
and supercritical regimes. Namely, to obtain the order of the mixing-time 
in Theorem 2 (critical window) , we study the magnetization chain using the 
arguments that appear in the proof of Theorem 1 (high temperature regime) . 
It is then straightforward to show that the mixing-time of the entire Glauber 
dynamics has the very same order. In turn, the spectral-gap in the critical 
window is obtained using arguments similar to those used in the proof of 
Theorem 3 (low temperature regime). In light of this, the following sketch 
will focus on the two non-critical temperature regimes. 

2.1. High temperature regime. 

Upper bound for mixing. As mentioned above, a key element in the proof is 
the analysis of the normalized magnetization chain, (St), which is the aver- 
age spin in the system. That is, for a given configuration a, we define S(a) 
to be i Yli a (^)i an d it is easy to verify that this function of the dynamics is 
an irreducible and aperiodic Markov chain. Clearly, a necessary condition 
for the mixing of the dynamics is the mixing of its magnetization, but inter- 
estingly, in our case the converse essentially holds as well. For instance, as 
we later explain, in the special case where the starting state is the all-plus 
configuration, by symmetry these two chains have precisely the same total 
variation distance from equilibrium at any given time. 

In order to determine the behavior of the chain (St), we first keep track 
of its expected value along the Glauber dynamics. To simplify the sketch 
of the argument, suppose that our starting configuration is somewhere near 
the all-plus configuration. In this case, one can show that ~ESt is monotone 
decreasing in t, and drops to order \J\jbn precisely at the cutoff point. 
Moreover, if we allow the dynamics to perform another 0(n/<5) steps (our 
cutoff window), then the magnetization will hit (or depending on the 
parity of n) with probability arbitrarily close to 1. At that point, we essen- 
tially achieve the mixing of the magnetization chain. 
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It remains to extend the mixing of the magnetization chain to the mix- 
ing of the entire Glauber dynamics. Roughly, keeping in mind the above 
comment on the symmetric case of the all-plus starting configuration, one 
can apply a similar argument to an arbitrary starting configuration a, by 
separately treating the set of spins which were initially positive and those 
which were initially negative. Indeed, it was shown in [11] that the following 
holds for j3 < 1 fixed (strictly subcritical regime). After a "burn-in" period 
of order n steps, the magnetization typically becomes not too biased. Next, 
if one runs two instances of the dynamics, from two such starting config- 
urations (where the magnetization is not too biased), then by the time it 
takes their magnetization chains to coalesce, the entire configurations be- 
come relatively similar. This was established by a so-called Two Coordinate 
Chain analysis, where the two coordinates correspond to the current sum 
of spins along the set of sites which were initially either positive or negative 
respectively. 

By extending the above Two Coordinate Chain Theorem to the case of 
(3 = 1 — 5 where 5 = 5(n) satisfies 5 2 n — > oo, and combining it with second 
moment arguments and some additional ideas, we were able to show that 
the above behavior holds throughout this mildly subcritical regime. The 
burn-in time required for the typical magnetization to become "balanced" 
now has order n/S, and so does the time it takes the full dynamics of two 
chains to coalesce once their magnetization chains have coalesced. Thus, 
these two periods are conveniently absorbed in our cutoff window, making 
the cutoff of the magnetization chain the dominant factor in the mixing of 
the entire Glauber dynamics. 

Lower bound for mixing. While the above mentioned Two Coordinate Chain 
analysis was required in order to show that the entire Glauber dynamics 
mixes fairly quickly once its magnetization chain reaches equilibrium, the 
converse is immediate. Thus, we will deduce the lower bound on the mixing 
time of the dynamics solely from its magnetization chain. 

The upper bound in this regime relied on an analysis of the first and 
second moments of the magnetization chain, however this approach is too 
coarse to provide a precise lower bound for the cutoff. We therefore resort 
to establishing an upper bound on the third moment of the magnetization 
chain, using which we are able to fine-tune our analysis of how its first mo- 
ment changes along time. Examining the state of the system order n/S 
steps before the alleged cutoff point, using concentration inequalities, we 
show that the magnetization chain is typically substantially far from 0. Re- 
calling Figure 2, this implies a lower bound on the total variation distance 
of the magnetization chain to stationarity, as required. 



8 



JIAN DING, EYAL LUBETZKY AND YUVAL PERES 



Spectral gap analysis. In the previous arguments, we stated that the magne- 
tization chain essentially dominates the mixing-time of the entire dynamics. 
An even stronger statement holds for the spectral gap: the Glauber dynam- 
ics and its magnetization chain have precisely the same spectral gap, and 
it is in both cases attained by the second largest eigenvalue. We therefore 
turn to establish the spectral gap of (St). 

The lower bound follows directly from the contraction properties of the 
chain in this regime. To obtain a matching upper bound, we use the Dirichlet 
representation for the spectral gap, combined with an appropriate bound on 
the fourth moment of the magnetization chain. 

2.2. Low temperature regime. 

Exponential mixing. As mentioned above, the exponential mixing in this 
regime follows directly from the behavior of the magnetization chain, which 
has a bottleneck between ±£. To show this, we analyze the effective resis- 
tance between these two centers of mass, and obtain the precise order of the 
commute time between them. Additional arguments show that the mixing 
time of the entire Glauber dynamics in this regime has the same order. 

Spectral gap analysis. In the above mentioned proof of the exponential mix- 
ing, we establish that the commute time of the magnetization chain between 
and C has the same order as the hitting time from 1 to 0. We can therefore 
apply a recent result of [5] for general birth- and-death chains, which implies 
that in this case the inverse of the spectral-gap (known as the relaxation- 
time) and the mixing-time must have the same order. 

3. Preliminaries 

3.1. The magnetization chain. The normalized magnetization of a con- 
figuration o~ E Q, denoted by S(a), is defined as 



Suppose that the current state of the Glauber dynamics is a, and that site i 
has been selected to have its spin updated. By definition, the probability of 
updating this site to a positive spin is given by p + (S(a) — a(i)/n), where 



Similarly, the probability of updating the spin of site i to a negative one is 
given by p~ (S(a) — a(i)/n), where 




p + (s) : 




(3.1) 



e f3s _|_ e -/3s 



2 




(3.2) 
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It follows that the (normalized) magnetization of the Glauber dynamics at 
each step is a Markov chain, (St), with the following transition kernel: 

ri±-p- (a _ n -i) xj = 8 -i 

Pm(s,s') = | i=*p+( a + n -i) if s' = S + |, 

(3.3) 

An immediate important property that the above reveals is the symmetry of 
St', the distribution of (St+i \ St = s) is precisely that of (—St+i \ St = —s). 

As evident from the above transition rules, the behavior of the Hyperbolic 
tangent will be useful in many arguments. This is illustrated in the following 
simple calculation, showing that the minimum over the holding probabili- 
ties of the magnetization chain is nearly |. Indeed, since the derivative of 
tanh(ar) is bounded away from and 1 for all x 6 [0,(3] and any j3 = 0(1), 
the Mean Value Theorem gives 

Pm(s,s+D = ^ (1 + tanhOfo)) + 0(n~ l ) , 
Pm(s,s-D = ^ (1 - tanhOfo)) + 0(n~ l ) , (3.4) 
Pm(s,s) = i (1 + stanh(/3s)) - Ofa" 1 ) . 

Therefore, the holding probability in state s is at least \ — O(^). In fact, 
since tanh(x) is monotone increasing, Pm(s,s) < \ + ^stanh(/3s) for all s, 
hence these probability are also bounded from above by |(1 + tanh(/3)) < 1. 

Using the above fact, the next lemma will provide an upper bound for 
the coalescence time of two magnetization chains, St and St, in terms of the 
hitting time To, defined as To := min{t : \St\ < n" 1 }. 

Lemma 3.1. Let (St) and (St) denote two magnetization chains, started 
from two arbitrary states. Then for any e > there exists some c £ > 0, 
such that the following holds: ifT > satisfies Pi(to > T) < e then St and 
St can be coupled in a way such that they coalesce within at most c £ T steps 
with probability at least 1 — e. 

Proof. Assume without loss of generality that \Sq\ < \Sq\, and by symmetry, 
that a = | S 1 > 0. Define 

r:=min{t: \S t \ < \S t \ + |} . 

Recalling the definition of To, clearly we must have t < tq. Next, since 
the holding probability of St at any state s is bounded away from and 1 
for large n (by the discussion preceding the lemma), there clearly exists a 
constant < b < 1 such that 

P (S t+1 = St+i | \s t - s t \ < l) > b > 



10 



JIAN DING, EYAL LUBETZKY AND YUVAL PERES 



(for instance, one may choose b = (1 — tanh(/3)) for a sufficiently large 
n). It therefore follows that |5V+i| = |5V+i| with probability at least b. 

Condition on this event. We claim that in this case, the coalescence of (St) 
and (St) (rather than just their absolute values) occurs at some t < tq + 1 
with probability at least b. The case SV+i = <SV+i is immediate, and it 
remains to deal with the case SV+i = — S T +i- Let us couple (St) and (St) 
so that the property St = —St is maintained henceforth. Thus, at time 
t = tq we obtain \St — St\ = 2\St\ < — , and with probability b this yields 
St+i = St+i- 

Clearly, our assumption on T and the fact that < a < 1 together give 

P CT (r > T) < P^ro >T)<e . 

Thus, with probability at least (1 — e)b 2 , the coalescence time of (St) and 
(St) is at most T. Repeating this experiment a sufficiently large number of 
times then completes the proof. ■ 

In order to establish cutoff for the magnetization chain (St), we will need 
to carefully track its moments along the Glauber dynamics. By definition 
(see (3.3)), the behavior of these moments is governed by the Hyperbolic 
tangent function, as demonstrated by the following useful form for the con- 
ditional expectation of St+i given St (see also [11, (2.13)]). 

E [St+i | S t = s] = (s + l)P M {s, s + l)+ sP M (s, s) + (s - l)P M (s, s - I) 
= (1 -rr l )s + ip(s) - iP(s) , (3.5) 

where 

<f(s) = ip(s, (3, n) := — [tanh ((3(s + n' 1 )) + tanh (/3(s - n" 1 ))] , 

iP(s) = ip(s, (3, n) := ^- [tanh (f3(s + n -1 )) - tanh ((3(s - n' 1 ))] . 

3.2. From magnetization equilibrium to full mixing. The motivation 
for studying the magnetization chain is that its mixing essentially dominates 
the full mixing of the Glauber dynamics. This is demonstrated by the next 
straightforward lemma (see also [11, Lemma 3.4]), which shows that in the 
special case where the starting point is the all-plus configuration, the mixing 
of the magnetization is precisely equivalent to that of the entire dynamics. 

Lemma 3.2. Let (Xt) be an instance of the Glauber dynamics for the mean 
field Ising model starting from the all-plus configuration, namely, o~q = 1, 
and let St = S(Xt) be its magnetization chain. Then 

\\Pl(X t €■)- M„||tV = ||Pl(S* G •) - TnllTV , (3-6) 
where 7r n is the stationary distribution of the magnetization chain. 
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Proof. For any s 6 {-1, -1 + £ , . . . , 1 - f , 1}, let U s := {a G U : 5(a) = s}. 
Since by symmetry, both fi n (- \ £l s ) and Pi(Af £ • | St = s) are uniformly 
distributed over Q s , the following holds: 



In the general case where the Glauber dynamics starts from an arbitrary 
configuration do, though the above equivalence (3.6) no longer holds, the 
magnetization still dominates the full mixing of the dynamics in the following 
sense. The full coalescence of two instances of the dynamics occurs within 
order n log n steps once the magnetization chains have coalesced. 

Lemma 3.3 ([11, Lemma 2.9]). Let a, a £ £1 be such that S(cr) = S(cr). For 
a coupling (At, Xt), define the coupling time t x x := min{i > : Xt = Xt}. 
Then for a sufficiently large cq > there exists a coupling (Xt,Xt) of the 
Glauber dynamics with initial states Xq = a and Xq = a such that 



Though Lemma 3.3 holds for any temperature, it will only prove useful in 
the critical and low temperature regimes. At high temperature, using more 
delicate arguments, we will establish full mixing within order of j steps once 
the magnetization chains have coalesced. That is, the extra steps required 
to achieve full mixing, once the magnetization chain cutoff had occurred, are 
absorbed in the cutoff window. Thus, in this regime, the entire dynamics has 
cutoff precisely when its magnetization chain does (with the same window). 

3.3. Contraction and one-dimensional Markov chains. We say that a 
Markov chain, assuming values in M, is contracting, if the expected distance 
between two chains after a single step decreases by some factor bounded 
away from 0. As we later show, the magnetization chain is contracting 
at high temperatures, a fact which will have several useful consequences. 
One example of this is the following straightforward lemma of [11], which 
provides a bound on the variance of the chain. Here and throughout the 
paper, the notation P z , E 2 and Var^ will denote the probability, expectation 
and variance respectively given that the starting state is z. 

Lemma 3.4 ([11, Lemma 2.6]). Let {Zt) be a Markov chain taking values 
in R and with transition matrix P. Suppose that there is some < p < 1 





Pi G •) - vr n || TV • 



lim sup P CT) 5- f t x £ > cqu log n j = . 
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such that for all pairs of starting states (z, z), 

\B z [Z t ]-B,[Z t ]\< p l \z-l\. (3.7) 

Then vt := sup zo Var^ (Z t ) satisfies vt < v\ min [t, 1/ (l — p 2 ) }. 

Remark. By following the original proof of the above lemma, one can readily 
extend it to the case p > 1 and get the following bound: 

vt<vi-p 2t min (M/Co 2 -!)} • (3.8) 

This bound will prove to be effective for reasonably small values of t in the 
critical window, where although the magnetization chain is not contracting, 
p is only slightly larger than 1. 

Another useful property of the magnetization chain in the high tempera- 
ture regime is its drift towards 0. As we later show, in this regime, for any 
s > we have E [St+i | = s] < s, and with probability bounded below by 
a constant we have St+i < St- We thus refer to the following lemma of [12]: 

Lemma 3.5 ([12, Chapter 18]). Let (Wi)t>o be a non-negative supermartin- 
gale and r be a stopping time such 

(i) W = k, 

(ii) W t+ i -W t <B, 

(hi) Var(W t+ i | Tt) > cr 2 > on the event r > t . 
Ifu> 4B 2 /(3a 2 ), then P fc (r > u) < 

This lemma, together with the above mentioned properties of (St), yields 
the following immediate corollary: 

Corollary 3.6 ([11, Lemma 2.5]). Let < 1, and suppose that n is even. 
There exists a constant c such that, for all s and for all u, t > 0, 

CTl S\ 

P( \S U \ > 0, . . . , \S u+t \ > I S u = s) < -jJ- . (3.9) 

Finally, our analysis of the spectral gap of the magnetization chain will 
require several results concerning birth-and-death chains from [5]. In what 
follows and throughout the paper, the relaxation-time of a chain, t REL , is 
defined to be gap"" 1 , where gap denotes its spectral-gap. We say that a 
chain is b-lazy if all its holding probabilities are at least b, or simply lazy for 
the useful case of b = ^. Finally, given an ergodic birth-and-death chain on 
X = {0, 1, . . . , n} with stationary distribution ir, the quantile state Q(a), for 
< a < 1, is defined to be the smallest i £ X such that 7r({0, . . . ,i}) > a. 
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Lemma 3.7 ([5, Lemma 2.9]). Let X(t) be a lazy irreducible birth-and- 
death chain on {0, 1, . . . , n}, and suppose that for some < e < we have 
£rel < £ 4 • EoTq( 1 _ £ - ) . Then for any fixed e<a<(3<l — e: 

Eq( q )T Q(/3) < — y^W • E 0T Q( 1) • (3.10) 

Lemma 3.8 ([5, Lemma 2.3]). For any fixed < e < 1 and lazy irreducible 
birth- and- death chain X, the following holds for any t: 

||P*(0, •) - vt||tv < Po(r Q( i_ £) > t) + e , (3.11) 

and for all k E O, 

||P*(&, •) - ttHtv < P fc (max{r Q(£) ,rQ (1 _ e) } > t) + 2e . (3.12) 

Remark. As argued in [5] (see Theorem 3.1 and its proof), the above two 
lemmas also hold for the case where the birth-and-death chain is not lazy 
but rather Mazy for some constant b > 0. The formulation for this more 
general case incurs a cost of a slightly different constant in (3.10), and re- 
placing t with t/C (for some constant C) in (3.11) and (3.12). As we already 
established (recall (3.4)), the magnetization chain is indeed Mazy for any 
constant b < \ and a sufficiently large n. 

3.4. Monotone coupling. A useful tool throughout our arguments is the 
monotone coupling of two instances of the Glauber dynamics (Xt) and (Xf), 
which maintains a coordinate-wise inequality between the corresponding 
configurations. That is, given two configurations a > a (i.e., a(i) > a(i) for 
all i), it is possible to generate the next two states a' and a' by updating the 
same site in both, in a manner that ensures that a' > a 1 . More precisely, 
we draw a random variable I uniformly over {1,2, ... ,n} and independently 
draw another random variable U uniformly over [0,1]. To generate a' from 
a, we update site / to +1 if U < p + (jS(a) — ^pjj otherwise a' (I) = — 1. 
We perform an analogous process in order to generate a' from a, using the 
same I and U as before. The monotonicity of the function p + guarantees 
that a' > a', and by repeating this process, we obtain a coupling of the two 
instances of the Glauber dynamics that always maintains monotonicity. 

Clearly, the above coupling induces a monotone coupling for the two cor- 
responding magnetization chains. We say that a birth-and-death chain with 
a transition kernel P and a state-space X = {0,1,..., n} is monotone if 
P(i,i + 1) + P(i + < 1 for every i < n. It is easy to verify that this 
condition is equivalent to the existence of a monotone coupling, and that for 
such a chain, if / : X — > K is a monotone increasing (decreasing) function 
then so is Pf (see, e.g., [5, Lemma 4.1]). 
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3.5. The spectral gap of the dynamics and its magnetization chain. 

To analyze the spectral gap of the Glauber dynamics, we establish the fol- 
lowing lemma which reduces this problem to determining the spectral-gap 
of the one-dimensional magnetization chain. Its proof relies on increasing 
eigenfunctions, following the ideas of [13]. 

Proposition 3.9. The Glauber dynamics for the mean-field Ising model 
and its one- dimensional magnetization chain have the same spectral gap. 
Furthermore, both gaps are attained by the largest nontrivial eigenvalue. 

Proof. We will first show that the one-dimensional magnetization chain has 
an increasing eigenfunction, corresponding to the second eigenvalue. 

Recalling that St assumes values in X := {—1, — 1 + j|, . . . , 1 — \, 1}, let 
M denote its transition matrix, and let ir denote its stationary distribution. 
Let I = Oq > 8i > . . . > 9 n be the n+1 eigenvalues of M, corresponding to 
the eigenfunctions /o = 1, /i, . . . , f n . Define 9 = max{#i, \0 n \\, and notice 
that, as St is aperiodic and irreducible, < 9 < 1. Furthermore, by the 
existence of the monotone coupling for St and the discussion in the previous 
subsection, whenever a function / : X — > IR is increasing so is Mf. 

Define / : / — > R by / := f\ + / n + Kl, where 1 is the identity function 
and K > is sufficiently large to ensure that / is monotone increasing 
(e.g., K = f ||/i + /n||i°° easily suffices). Notice that, by symmetry of St, 
7t(x) = 7r(— x) for all x £ X, and in particular YlxeX X7T ( X ) = 0, that is to 
say, (1, / ) i2(7r) = 0. Recalling that for all i ^ j we have (/j, fj) L 2^ = 0, it 
follows that for some qi, . . . ,q n G K we have / = Y17=l w *th ?i / and 
q n 0, and thus 

n 

(e-i M ) m f = Y^q t (9 l /9rf i . 
i=l 

Next, define 

f qifl if 6 = 0! , , / q n fn H9 = -9 n 

g = < , and h = < , 

[ otherwise [ otherwise 

and notice that 

lim (9~ 1 M) 2m f = g + h , and lim (9' 1 M) 2m+1 f = g - h . 

As stated above, M m f is increasing for all m, and thus so are the two limits 
g + h and g—h above, as well as their sum. We deduce that g is an increasing 
function, and next claim that g ^ 0. Indeed, if g = then both h and —h 
are increasing functions, hence necessarily h = as well; this would imply 
that q± = q n = 0, thus contradicting our construction of /. 
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We deduce that g is an increasing eigenfunction corresponding to 9\ = 9, 
and next wish to show that it is strictly increasing. Recall that for all x E X, 

(Mg)(x) = M(x,x-^jg(x--^j +M(x, x)g{x) + m(x,x+ -^jg(x+ . 

Therefore, if for some x S X we had g(x — §) = s( x ) — 0, the fact that g is 
increasing would imply that 

9 ig (x) = (Mg)(x) > g(x) > , 

and analogously, if g(x) = g(x + |) < we could write 

9 ig (x) = (Mg)(x) < g(x) < . 

In either case, since < 9\ < 1 (recall that 9\ = 0) this would in turn lead 
to g{x) = 0. By inductively substituting this fact in the above equation for 
(Mg)(x), we would immediately get g = 0, a contradiction. 

Let 1 = Ao > Ai > . . . > A^i^ denote the eigenvalues of the Glauber 
dynamics, and let A := max{Ai, |A2"-i|}. We translate g into a function 
G : Q — > M in the obvious manner: 

1 n 

G{a) :=g(S(a))=g[-J2^) ■ 

i=i 

One can verify that G is indeed an eigenfunction of the Glauber dynamics 
corresponding to the eigenvalue 9i, and clearly G is strictly increasing with 
respect to the coordinate- wise partial order on Q. At this point, we refer to 
the following lemma of [13]: 

Lemma 3.10 ([13, Lemma 4]). Let P be the transition matrix of the Glauber 
dynamics, and let Xi be its second largest eigenvalue. If P has a strictly 
increasing eigenfunction f , then f corresponds to X\. 

The above lemma immediately implies that G corresponds to the second 
eigenvalue of Glauber dynamics, which we denote by Ax, and thus Ai = 9\. 

It remains to show that A = Ai. To see this, first recall that all the holding 
probabilities of St are bounded away from 0, and the same applies to the 
entire Glauber dynamics by definition (the magnetization remains the same 
if and only if the configuration remains the same). Therefore, both 9 n and 
A2«-i are bounded away from —1, and it remains to show that gap = o(l) 
for the Glauber dynamics (and hence also for its magnetization chain). 

To see this, suppose P is the transition kernel of the Glauber dynam- 
ics, and recall the Dirichlet representation for the second eigenvalue of a 
reversible chain (see [12, Lemma 13.7], and also [1, Chapter 3]): 

1-Ai = min{^^- : /#0, E^(/) = o} , (3.13) 
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where E Mn (/) denotes and 



S(f) = ((/ -P)f,f)^ = - £ [/(a) - /(a')] 2 /i»P(<x, a') . 



By considering the sum of spins, h(a) = Y17=l a (^)' we ^(^) — 2, and 
since the spins are positively correlated, Var Aln £V <r(i) > n. It follows that 



1 - Ai < 2/n , 

and thus gap = 1 — Ai = 1 — 6\ for both the Glauber dynamics and its 



In this section we prove Theorem 1. Subsection 4.1 establishes the cutoff 
of the magnetization chain, which immediately provides a lower bound on 
the mixing time of the entire dynamics. The matching upper bound, which 
completes the proof of cutoff for the Glauber dynamics, is given in Subsection 
4.2. The spectral gap analysis appears in Subsection 4.3. Unless stated 
otherwise, assume throughout this section that (5 = 1 — 5 where 5 2 n — ► oo. 

4.1. Cutoff for the magnetization chain. Clearly, the mixing of the 
Glauber dynamics ensures the mixing of its magnetization. Interestingly, 
the converse is also essentially true, as the mixing of the magnetization 
turns out to be the most significant part in the mixing of the full Glauber 
dynamics. We thus wish to prove the following cutoff result: 

Theorem 4.1. Let (3 = 1 — 5, where 5 > satisfies 5 2 n — ► oo. Then the 
corresponding magnetization chain (St) exhibits cutoff at time \ ■ j log(<5 2 n) 
with a window of order n/5. 

Notice that Lemma 3.2 then gives the following corollary for the special 
case where the initial state of the dynamics is the all-plus configuration: 

Corollary 4.2. Let 5 = 5(n) > be such that 5 2 n — ► oo with n, and let (Xt) 
denote the Glauber dynamics for the mean-field Ising model with parameter 
(3 = 1 — 5, started from the all-plus configuration. Then (Xt) exhibits cutoff 
at time ^(n/5) log(5 2 n) with window size n/5. 

4.1.1. Upper bound. Our goal in this subsection is to show the following: 



where d n (-) is with respect to the magnetization chain (St) and its stationary 
distribution. This will be obtained using an upper bound on the coalescence 
time of two instances of the magnetization chain. Given the properties of 



magnetization chain, as required. 



4. High temperature regime 




(4.1) 
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its stationary distribution (see Figure 2), we will mainly be interested in the 
time it takes this chain to hit near 0. The following theorem provides an 
upper bound for that hitting time. 

Theorem 4.3. For < (3 < l + 0(n~ 1 ^ 2 ), consider the magnetization chain 
started from some arbitrary state sq, and let tq = min{t : \St\ < n -1 }. Write 
(3 = 1 — 5, and for 7 > define 



log(5 2 n) + (7 + 3) - 5 2 n -► 00 , 

26 5 (42) 



t n (j) 

'2OO + 67 (l + 6v / S 2 n)) n 3/2 5 2 n = 0{l) . 

Then there exists some c > such that P So (to > t n (7)) < c l \/l ■ 

Proof. For any t > 1, define: 

st ■= E so [|St|l{ r()>t }] . 

Suppose s > 0. Recalling (3.5) and bearing in mind the concavity of the 
Hyperbolic tangent and the fact that f/>(s) > 0, we obtain that 

E(5t+i I S t = s) < s + - ( tanh(/3s) - s) . 
Using symmetry for the case s < 0, we can then deduce that 

E [\S t+ i\ I St] < \S t \ + i(tanh(/9|5t|) - \S t \) for any t < r . (4.3) 

Hence, combining the concavity of the Hyperbolic tangent together with 
Jensen's inequality yields 

s t +i < (l--)st + - tanh(/3s 4 ) . (4.4) 
\ n J n 

Since the Taylor expansion of tanh(x) is 

x^ 23?^ 17 ' x^ 
tanh(x) = x ~Y + ~l5~Yl5 + ' (4 " 5) 

we have tanh(x) < x — ^~ for < x < 1, giving 



st+i < ( 1 - -) s t + - tanh(/3s f ) < (l - -) s t + -f3s t 
\ nJ n \ nJ n 



n) " ' n v "' V nJ " ' n 5n 

= «-^-%£. (4.6) 
n an 

For some 1 < a < 2 to be defined later, set 

b{ = a~~ % /4 , and Uj = min{t : st < 6j} . 
Notice that st is decreasing in t by (4.6), thus for every t G [uj, lii+i] we have 

6j/o = 6 i+1 < s t < 6j . 
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It follows that 

5 h bf 

s t +i < s t , 

n a ocrra 

and 

fa-l\,(5bi bf \ 5(a-l)a 2 n w . 

For the case <5 2 n — ► 00, define: 

io = min{i : 6j < 1/ vJra} . 

The following holds: 

E'° . . 5(a — l)a 2 n s-^5(a — l)a 2 n v-^ 5(o — l)a 2 ra 

<^- , ^5:-b?Tir*E- L y— +E-4*?— 

6 2 >5 6? < 5 

5na 2 a — 1 n , , „ 9 . 

^ AmiTT + ^ 7 log <5 2 n , 

o(a + 1) 2 log a 

where in the last inequality we used the fact that the series {b^ 2 } is a 
geometric series with a ratio a 2 , and that, as bf > l/(6n) for all i < the 
number of summands such that bf < 5 is at most log a (V <5 2 n). Therefore, 
choosing a = 1 + n _1 , we deduce that: 

X> + i - m) < (J + Otn- 1 )) 2 + Q + {n^ - 6 log(5 2 n) 

<3™ + ^log(5 2 n) , (4.8) 

where the last inequality holds for any sufficiently large n. Combining the 
above inequality and the definition of zq, we deduce that 

( s t „ (0 ) = ) E S0 [|S tn(0) |l {ro>tn(0)} ] < . (4.9) 

Thus, by Corollary 3.6 (after taking expectation), for some fixed c > 

P(r > t„(7)) < • 
For the case 5 2 n = O(l), choose a = 2, that is, b{ = 2~(* +2 ), and define 
h = min{i : bi < n* 1 / 4 V 5y/\6\} . (4.10) 
Substituting a = 2 in (4.7), while noting that 5 > —^bf for all i < i±, gives 

Ev — v 2un \ — v n n 
(u i+1 - «o < £ ^TV2 < 100 E ^ 20 V 

j=l i=l * i=l 1 'l 

n 



< ( 200n 3/2 A 8— ] < 20()» :! 2 
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in 



where the last inequality in the first line incorporated the geometric sum 
over {br 2 }. By (4.10), 

b h < n- 1 / 4 V 5^/\6\ < n- l ' A (l + ^n) 1 / 4 ^ , 

and as in the subcritical case, we now combine the above results with an 
application of Corollary 3.6 (after taking expectation), and deduce that for 
some absolute constant c > 0, 



P(r > t«(7)) < , 



as required. 



Apart from drifting toward 0, and as we had previously mentioned, the 
magnetization chain at high temperatures is in fact contracting; this is a 
special case of the following lemma. 

Lemma 4.4. Let (St) and (St) be the corresponding magnetization chains 
of two instances of the Glauber dynamics for some (3 = 1 — 8 ( where 5 is not 
necessarily positive), and put Dt := St — St- The following then holds: 



E[A+i - D t | D t ] < — D t + 



\D t 



+ 0(n 



n n 

Proof. By definition (recall (3.5)), we have 

E[A+i - D t | D t ] = HSt+i -St + St- St+i I D t ] 

St -s t 



(4.11) 



+ 



ii 



tp(St) - <p(St) - ^St)-^(St) 



The Mean Value Theorem implies that 



ip(St) - p(St) < -(St - s t ) 



n 



and applying Taylor expansions on tanh(x) around (5 St and (3 St, we deduce 
that 



V>(s*)-v(St) 



St 



n 2 cosh 2 (f3S t ) n 2 cosh 2 ((3 S t ) 



+ o 



1 



ii' 



Since the derivative of the function xj cosh 2 (/3x) is bounded by 1, another 
application of the Mean Value Theorem gives 



r/»(S t )-V(S t ) 
Altogether, we obtain (4.11), as required. 



< lS ^ + o(\ 
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Indeed, the above lemma ensures that in the high temperature regime, 
(3 = 1 — 5 where 5 > 0, the magnetization chain is contracting: 

E [\D t+1 \ I D t ) <(l ) \D t \ for any sufficiently large n . (4.12) 



2n, 

We are now ready to prove that hitting near essentially ensures the mixing 
of the magnetization. 

Lemma 4.5. Let (3 = 1 — 5 for 5 > with 5 2 n — > oo, (X t ) and (X t ) 
be two instances of the dynamics started from arbitrary states o~q and 5"o 
respectively, and (St) and (St) be their corresponding magnetization chains. 
Let r mag denote the coalescence time r mag := min{i : St = St}, and t n (^) be 
as defined in Theorem 4-3. Then there exists some constant c > such that 

P (r mag > t„(3 7 )) < c/V7 for all 7 > . (4.13) 

Proof. Set T = i n ( 7 )- We claim that the following holds for large n: 

\ES t \ < -== and |E5 t | < -j= for all i > T . (4.14) 
V Sn V on 



To see this, first consider the case where n is even. The above inequality 
then follows directly from (4.9) and the decreasing property of St (see (4.6)), 
combined with the fact that EoS* = (and thus ESt = for all t > to). In 
fact, in case n is even, |ESj| and |ESj| are both at most \ j\fbn for all t > T. 
For the case of n odd (where there is no state for the magnetization chain, 
and To is the hitting time to ±^), a simple way to show that (4.14) holds is to 
bound |Ei5t|. By definition, P M (\, \) > P M (L -£) (see (3.3)). Combined 

n 

with the symmetry of the positive and negative parts of the magnetization 
chain, one can then verify by induction that -P/t/^n'n) — ^AfCn' - n) ^ or 
any odd k > and any t. Therefore, by symmetry as well as the fact that 
E S0 S£ < so f° r positive so, we conclude that |Ei5j| is decreasing with t, and 

n 

thus is bounded by — . This implies that (4.14) holds for odd n as well. 

Combining (4.14) with the Cauchy-Schwartz inequality we obtain that for 
any t >T 



E\S t - St\ < E\St\ + E\St\ < d Var(5 4 ) + i- + J V&i(S t ) + ^ . 

Now, combining Lemma 3.4 and Lemma 4.4 (and in particular, (4.12)), we 
deduce that Var St < tz:, and plugging this into the above inequality gives 

E\S t - S t \ < 4^= for any t > T . 
V 5n 



We next wish to show that within 2^n/5 additional steps, St and St coalesce 
with probability at least 1 — for some constant c > 0. 
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Consider time T, and let D t := St — St- Recall that we have already 
established that 

ED T < lO/Vfa , (4.15) 

and assume without loss of generality that Dt > 0. We now run the mag- 
netization chains St and St independently for T < t < t\, where 

T i: =min{t>T:AG{0,-f}} , 

and let Ft be the cr-field generated by these two chains up to time t. By 
Lemma 4.4, we deduce that for sufficiently large values of n, if Dt > then 

E[A+i-A|ff]<- r A<0, (4.16) 

In 

and Dt is a supermartingale with respect to T%. Hence, so is 

n 

W t := D T+t ■ -l {T1>t } , 

and it is easy to verify that Wt satisfies the conditions of Lemma 3.5 (recall 
the upper bound on the holding probability of the magnetization chain, 
as well as the fact that at most one spin is updated at any given step). 
Therefore, for some constant c > 0, 

P (n > t n (2 7 ) | D T ) = P(W > 0, W 1 > 0, . . . , W tn{ ^)- T > \ D T ) 
^ cuDt 
~ 's/yn/S 

Taking expectation and plugging in (4.15), we get that for some constant c', 

P (n > t n (2 7 )) < — . (4.17) 
V7 



From time t\ and onward, we couple St and St using a monotone coupling, 
thus Dt becomes a non-negative supermartingale with D n < ^. By (4.16), 

E [A+i - D t | F t ] < . for n < t < r mag , 

and therefore, the Optional Stopping Theorem for non-negative supermartin- 
gales implies that, for some constant c", 

P (r mag - n > n/S) < E(Tmas / 7 Tl) < - . (4.18) 

n/o 7 

Combining (4.17) and (4.18) we deduce that for some constant c, 

P (T mag > ^ (37)) < — , 

completing the proof. ■ 
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4.1.2. Lower bound. We need to prove that the following statement holds 
for the distance of the magnetization at time t from stationarity: 

- • — log(5 2 n) -7-1=1. (4.19) 
2 J 

The idea is to show that, at time g ■ 7 log(<5 2 n) — 7^, the expected magne- 
tization remains large. Standard concentration inequalities will then imply 
that the magnetization will typically be significantly far from 0, unlike its 
stationary distribution. 

To this end, we shall first analyze the third moment of the magnetization 
chain. Recalling the transition rule (3.3) of St under the notations (3.1), (3. 2) 

, 1 + tanh(/?s) , . 1 — tanh(/?s) 
P (*) = j ' P W = 2 ' 

the following holds: 
E [Sf +1 \S t = s] 

;) 3 + ^ + (, + n-)( s+ £) 

+ (l - l ~^P~( s - n ~ X ) - ^P + ( s + s 3 
6s 2 1 / 

= s 3 H [ - 2s + tanh (/3(s - n -1 )) + tanh (j3{s + n' 1 )) 

+ s (tanh ((3(s - n' 1 )) - tanh (/3(s + n -1 )) ) + ci^ + -^f . (4.20) 

As tanh(x) < x for x > 0, for every s > we get 

3 s 2 s co 

E [S 3 +1 I ft = s] < s 3 + — (-2s + /?(, - n" 1 ) + /?(* + n" 1 )) + c x ~ 2 + J 

= s 3_ 3 £ s 3 + CI + C2 (421) 

n n z n s 

If s = 0, the above also holds, since in that case |St_|_i| 3 < (2/n) 3 . Finally, 
by symmetry, if s < then the distribution of |5' 3 +1 | = — Sf +1 given St = s 
is the same as that of Sf +1 given St = \s\, and altogether we get: 

E[|5 m | 3 |ft = S ]<| S | 3 -3-| S | 3 + ^| S | + ^|. 
We deduce that 

E|5 t+1 | 3 <Ef|ft| 3 -3-|5 t | 3 + ^|ft| + ^ 



<(l--")E|ft| 3 + ^E|5 t | + ^. (4.22) 
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Note that the following statement holds for the first moment of St'. 
E so [\S t \] < ^(E S0 S t ) 2 + Var S0 (S*) 



so 16 ( S\* . , 4 
< A K) 2 + t-< 1-- M + 



Hence, 

E, |S m | 3 < ( 1 - -) E S0 |5 t | 3 + ^ f 1 - *Y N + ^= + ^ 
V n / n 2 \ nj n 2 \/Sn n d 

^ t ci, 4 c 2 <5 2 



^ol^t| T'/ 2 |^U| "I" 2 A? ^ „2 ' 

n n A y on n 

where rj = 1 — 5/n, and the extra error term involving c' 2 absorbs the change 

J J X2 \ * 



of coefficient of E So |Sj| 3 and also the 1/n 3 term. Iterating, we obtain 



EjSt+1 |'<^| S0 |3 + ^ WE ^ + (^ + ^) E 

j=0 \u Wit, / . =Q 

,3ti„ |3 , J C 1 M , ( A , c 2<^ 1 



ra 2 1 — rj 2 \n 2 \f5n n 2 J 1 — rf" 

< V 3t \so\ 3 + rf~ \^ + tAt2 + —- ( 423 ) 
on [on) 6 ^ n 

Define Z t := |5t|r/ - *, whence Zo = |5o| = |so] - Recalling (3.5), and com- 
bining the Taylor expansion of tanh(x) given in (4.5) with the fact that 
\ip(s)\ = O (s/n 2 ), we get that for s > 

3 

E[|5 m ||5 t = S ] >^-i__4 . 
L 1 J 2n n z 

By symmetry, an analogous statement holds for s < 0, and altogether we 
obtain that 

V[\S t+1 \\S t ]>n\S t \-^-^ . (4.24) 

Remark. Note that (4.24) in fact holds for any temperature, having followed 
from the basic definition of the transition rule of (St), rather than from any 
special properties that this chain may have in the high temperature regime. 

Rearranging the terms and multiplying by rf~ we obtain that for any 
sufficiently large n, 



E 



2 



1 



l-^)Z t - Zt+i | S t < -V^lSti 
n A J J n 
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where we used the fact that rj^ 1 < 2. Taking expectation and plugging in 
(4.23), we deduce that 



1 n ) Zt — Zf+3 



n \ on n 



(4.25) 
Set 

* = 25 log ^ 2n ) ~ in / 5 ' 
and notice that when n is sufficiently large, (l — \) ( t+is > < 2 for any t <t. 
Therefore, multiplying (4.25) by (1 — ^-)~*- t+1 ^ and summing over gives: 

2|s | 3 - ci 2r/-* / ci c' 2 5 



\s \ 2E S0 Z- t < n{l _ r]2) +t 5n 2\ s *\ + n{l _ rj) {( Sn )3/2 + n 

< 2|s | 3 + ci log(£ 2 re) + 4 + c 2 + 45 
~~ 5 25 2 n S 3 / 2 n 5rfi/ 2 ^fn 

= 2 \^l + o(V~5 + \s \), 

where the last inequality follows from the assumption 5 2 n — ► 00. We now 
select so = VS/3, which gives 

Vd/3 - 2E SQ Zj < 2v^/27 + o(V5) , 

and for a sufficiently large n we get 

E S0 Zj> V5/9 . 

Recalling the definition of Zt, and using the well known fact that (1 — x) > 
exp(— x/(l — x)) for < x < 1, we get that for a sufficiently large n, 

E S0 |Sf| > rfVS/9 > =: L . (4.26) 

lOv on 



Lemma 3.4 implies that max{Var So (S'i), Var /in (S't)} < 16/ Sn. Therefore, 
recalling that E Mn Sy = 0, Chebyshev's inequality gives 

16/(oV) _ 7 



P S0 (|5 f | < L/2) < P S0 (||5 f | - E S0 |5 ? || > L/2) < 
P Mn (|^| > L/2) < 



L 2 /4 

16/(oV) _ 



L 2 /4 ' 

Hence, letting 7r denote the stationary distribution of St, and ^ denote the 
set [— j, 4] , we obtain that 

||P S0 (5 f e •) - 7t||tv > vr(,4 L ) - P So (|5 f | G A x ) > 1 - 2ce^ , 

which immediately gives (4.19). ■ 
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4.2. Full Mixing of the Glauber dynamics. In order to boost the mix- 
ing of the magnetization into the full mixing of the configurations, we will 
need the following result, which was implicitly proved in [11, Sections 3.3, 
3.4] using a Two Coordinate Chain analysis. Although the authors of [11] 
were considering the case of < (3 < 1 fixed, one can follow the same argu- 
ments and extend this result to any (3 < 1. Following is this generalization 
of their result: 

Theorem 4.6 ([11]). Let (Xt) be an instance of the Glauber dynamics and 
\x n the stationary distribution of the dynamics. Suppose Xq is supported by 

n := {a G n : \S(a)\ < 1/2} . 

For any o~q G £lo an d & £ £1, we consider the dynamics (Xt) starting from 
o"o and an additional Glauber dynamics (Xt) starting from a, and define: 

r mag := min{t : S(X t ) = S(X t )} , 

U(a) := |{t : a(i) = a (i) = 1}\ , V(a) := \{i : a(i) = a (i) = -1}\ , 
3 := {a : mm{U(a), U(a ) - U(a), V(a),V(a ) - V(a))} > n/20} , 

R(t) := \u(X t )-U(X t )\ , 

Hi(t) := {r mag < t} , H 2 (h,t 2 ) := n% h {Xi G 3 A X { G E} . 
Then for any possible coupling of Xt and Xt, the following holds: 



max ||P CTo (X r2 G •) - ^nlkv < max P^^H^ri)) 

sen 



+ P ao A R n > cty/n/6) + PaoA H 2(n,r 2 )) + 



ac± n 



V r 2 - n V 8 



(4.27) 



where n < r 2 and a > 0. 



The rest of this subsection will be devoted to establishing a series of 
properties satisfied by the magnetization throughout the mildly subcritical 
case, in order to ultimately apply the above theorem. 

First, we shall show that any instance of the Glauber dynamics concen- 
trates on once it performs an initial burn-in period of n/5 steps. It 
suffices to show this for the dynamics started from so = 1: to see this, con- 
sider a monotone-coupling of the dynamics (Xt) starting from an arbitrary 
configuration, together with two additional instances of the dynamics, (X^~) 
starting from sq = 1 (from above) and (Xj~) starting from sq = — 1 (from 
below). By definition of the monotone-coupling, the chains (X^~) and (Xf) 
"trap" the chain (Xt), and by symmetry it indeed remains to show that 

Pl(|Sto| < 1/2) = 1 - o(l) , where to = n/5 . 
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Recalling (4.4), we have st+\ < (1 — ^)s< where st = E [|St|lr T0>t \] , thus 

Ei [\S to \l {ro>to} ] ^e" 1 . 

Adding this to the fact that EiSt l{To < to} = 0, which follows immediately 
from symmetry, we conclude that Ei5j < e _1 . Next, applying Lemma 3.4 
to our case and noting that (3.7) holds for p' = 1 — 4 (l — ntanh(^)) < 1 — 
we conclude that 

Var(S t ) < v£ < (-Y ? = ^ for all t . 
o \n J o on 

Hence, Chebyshev's inequality gives that \St \ < 1/2 with high probability. 
We may therefore assume henceforth that our initial configuration already 
belongs to some good state o"o £ fio- 
Next, set: 

T := t n (7) ,r := t n (2>y) ,n := ^(37) ,r 2 := t n (47) . 

We will next bound the terms in the righthand side of (4.27) in order. First, 
recall that Lemma 4.5 already provided us with a bound on the probability 
of H\{ri), by stating there for constant c > 

P(r mag > n) < . (4.28) 

Our next task is to provide an upper bound on R ri , and namely, to show 
that it typically has order at most \Jnf8. In order to obtain such a bound, 
we will analyze the sum of the spins over the set B := {i : (Tq{%) = 1}. Define 

ieB 

and consider the monotone-coupling of (Xt) with the chains (X^~) and (Xf) 
starting from the all-plus and all-minus positions respectively, such that 
X^ < Xt < Xf ' . By defining M t + and M t ~ accordingly, we get that 

E(M t (B)) 2 < E(M t + (S)) 2 + E{Mf{B)f = 2E(M+(5)) 2 . 

By (4.14), we immediately get that for t > T, \EMf(B)\ < y/f . We will 
next bound the variance of M^(B), by considering the following two cases: 

(i) If every pair of spins of X^ is positively correlated (since Xq is the all- 
plus configuration, by symmetry, the covariances of each pair of spins 
is the same), then we can infer that 

Var(M+(B)) < Var (1 £ X+(i)) = ^ Var (S(X+)) < ^ . 

ie[n] 
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(ii) Otherwise, every pair of spins of is negatively correlated, and it 
follows that 

Var(M+(S)) < ]T Var (^ + (*)) < 
Altogether, we conclude that for all t >T, 



n 
4 



E|M t (B)| < y E (M t (B)) < J 2 Var(M t + (_B) + 2 (EM t (B)) 



8n 2n n , d . 

£ Vt + t £8 vV (429) 

This immediately implies that 

~ ~ / ft 

ER ri = E\M n (B)- M ri (B)\ < E\M ri (B)\ + E\M ri (B)\ < 16 J - , 

and an application of Markov's inequality now gives 

P(i?n > a./?) < - • (4.30) 
V o a 



It remains to bound the probability of H2{r\, r<z). Define: 
Y:= Yl l{|M t (S)|>n/64}, 

ri<t<T2 

and notice that 

^ t=T\ ' 

Recall that the second inequality of (4.29) actually gives E\M t (B)\ 2 < ^f. 
Hence, a standard second moment argument gives 

1 



P{\M t {B)\ >n/64) = 0, 

on 



Altogether, F, ao Y = 0{1/ 5 2 ) and 

rJ\J{\M t (B)\>n/32A=o(^A . 

Applying an analogous argument to the chain {Xt), we obtain that 

rJ{j{mB)\>n/32})= (^ i ) , 

\ i =ri / \ / 

and combining the last two inequalities, we conclude that 

1 



P tT0 ,a[H 2 (r 1 ,r 2 ))=O[ ¥ ^) . (4.31) 
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Finally, we have established all the properties needed in order to apply 
Theorem 4.6. At the cost of a negligible number of burn-in steps, the state 
of (Xt) with high probability belongs to Oo- We may thus plug in (4.28), 
(4.30) and (4.31) into Theorem 4.6, choosing a = to obtain (4.1). 



4.3. Spectral gap Analysis. By Proposition 3.9, it suffices to determine 
the spectral gap of the magnetization chain. The lower bound will follow 
from the next lemma of [4] (see also [12, Theorem 13.1]) along with the 
contraction properties of the magnetization chain. 

Lemma 4.7 ([4]). Suppose 0, is a metric space with distance p. Let P be 
a transition matrix for a Markov chain, not necessarily reversible. Suppose 
there exists a constant 6 < 1 and for each x,y £ £1, there is a coupling 
(X^Yi) ofP(x,-) andP(y,-) satisfying 

E^MXi,*!)) < 9p{x,y) . 

If X is an eigenvalue of P different from 1, then |A| < 0. In particular, the 
spectral gap satisfies gap > 1 — 6. 

Recalling (4.12), the monotone coupling of St and St implies that 

r r 

E s , s |5i-5i| < (l- - + \s-s\ . 

Therefore, Lemma 4.7 ensures that gap > (1 + o(l))— . 

It remains to show a matching upper bound on gap, the spectral gap of 
the magnetization chain. Let M be the transition kernel of this chain, and tt 
be its stationary distribution. Similar to our final argument in Proposition 
3.9 (recall (3.13)), we apply the Dirichlet representation for the spectral gap 
(as given in [12, Lemma 13.7]) with respect to the function / being the 
identity map 1 on the space of normalized magnetization, we obtain that 

Eap£ (oyr (i.D. =1 — ^! — ' (4 ' 32) 

where E^S^ is the k-th moment of the stationary magnetization chain (St). 
Recall (4.24) (where rj = 1 — ^), and notice that the following slightly 
stronger inequality in fact holds: 

\St\ 3 \St\ 



E[sign(S t )St +1 \St} >ri\S, 



2n ri 



2 



(to see this, one needs to apply the same argument that led to (4.24), then 
verify the special cases S t S {0, i}). It thus follows that 

E [S t S t+1 | St] >r,S 2 t -^-%, 
L 1 J In n A 
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and plugging the above into (4.32) we get 

gap < --!--■ ^+^7. (1.*$) 



6_ 1_ B n Sf J_ 
n + 2~n ' + ^ ■ 



In order to bound the second term in (4.33), we need to give an upper 
bound for the fourth moment in terms of the second moment. The next 
argument is similar to the one used earlier to bound the third moment of 
the magnetization chain (see (4.20)), and hence will be described in a more 
concise manner. 

For convenience, we use the abbreviations h + := tanh m(s + n -1 )) and 
h~ := tanh (f3(s — n -1 )). By definition (see (3.3)) the following then holds: 

E[S t 4 +1 \S t = s] = s 4 + ^s 3 (-2s + h~ +h + + shT - h + ) 
+ ^s 2 (2 + h+ - h~ - sbT + h + ) 



„-s"(-2s + h +h + + sh - h + ) 



4 

+ —r (2 + h + - h~ - sh~ + h + ) 

( 45\ 4 12 2 16 
< 1 _ _ s 4 + 2 + 

\ n J n z n 4 

Now, taking expectation and letting the St be distributed according to tt, 
we obtain that 

Recalling that, as the spins are positively correlated, Var 7r (<S't) > — , we get 

e ^H)¥- (4 - 34) 

Plugging (4.34) into (4.33), we conclude that 

^( 1 + (i)) = (1+0(1) ^ 



5. The critical window 

In this section we prove Theorem 2, which establishes that the critical win- 
dow has a mixing-time of order n 3//2 without a cutoff, as well as a spectral- 
gap of order n~ 3 / 2 . 
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5.1. Upper bound. Let (Xt) denote the Glauber dynamics, started from 
an arbitrary configuration a, and let (Xt) denote the dynamics started from 
the stationary distribution \i n . As usual, let (St) and (St) denote the (nor- 
malized) magnetization chains of (Xt) and (Xt) respectively. 

Let e > 0. The case 5 2 n = 0(1) of Theorem 4.3 implies that, for a 
sufficiently large 7 > 0, P<j (tq > 7?i 3 / 2 ) < e. Plugging this into Lemma 
3.1, we deduce that there exists some c e > 0, such that the chains St and St 
coalesce after at most c £ n 3 ^ 2 steps with probability at least 1 — e. 

At this point, Lemma 3.3 implies that (Xt) and (Xt) coalesce after at most 
0(n 3 / 2 ) + 0(n log n) = 0(n 3 / 2 ) additional steps with probability arbitrarily 
close to 1, as required. 



5.2. Lower bound. Throughout this argument, recall that 5 is possibly 
negative, yet satisfies 5 2 n = 0(1). By (4.22), 

E\S t+1 f<v(\St\ 3 -3 6 -\Stf + ^\St\ + ^ 



n J n 2 n 3 

Recalling Lemma 4.4, and plugging the fact that 5 = 0(n -1 / 2 ) in (4.11), the 
following holds. If St and St are the magnetization chains corresponding to 
two instances of the Glauber dynamics, then for some constant c > and 
any sufficiently large n, 

E SjS -|^i - §i\ < (1 + cn- 3/2 )|s - s\ . (5.1) 

Combining this with the extended form of Lemma 3.4, as given in (3.8), we 
deduce that if t < en 3 / 2 for some small fixed e > 0, then Var So St < 4:t/n 2 . 
Therefore, 



E so [\S t \] < v/|E S0 5 t | 2 +Var S0 ^ < - ^ |s | + 



■n 



Therefore, 



E S0 |5 m r < I 1 - ^ E SoN 3 + °\ ( 1 - \s \ + ^ 
< v *E So \St\ s + i 1 t ^\so\ + ^ , 
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where again r) = 1 — S/n. Iterating, we obtain 

EjS m | 3 < v 3t \sof + X> 2j + ^ I> 3j 

i=o i=o 

o.. .0 . ci 77 2f_1 — 1 , , c\\/i r? 3 * — 1 
n z rj z — l n 6 T] 6 — 1 

< ?? 3t |so| 3 + ^4 - s 1 + ^V- -3i , (5.2) 

where the last inequality holds for sufficiently large n and t < en 3 / 2 with 
e > small enough (such a choice ensures that rj will be suitably small). 
Define Z t := |Si|?7 - *, whence Zq = \Sq\ = \sq\. Applying (4.24) (recall that 
it holds for any temperature) and using the fact that r/ _1 < 2, we get 

E[Z t+1 \S t )>Z t -^ (?T'N 3 + 0(l/n)) , 
for n large enough, hence 

E[Z t - Z t+1 \S t }<- {v'W + 0(l/n)) . 

n 

Taking expectation and plugging in (4.23), 

E S0 [Z t - Z t+1 ] < - ( V 2t \s \ 3 + °^\s \ + rf td ^- + 0{l/n) ) . (5.3) 

Set t = n 3//2 /A A for some large constant A such that ^ < if < 2. Summing 
over (5.3) we obtain that 

|s | - E S0 Z- t < )~ ^ | so | 3 + t 2 ^|s | + 1rf l • t b/2 /n A + 0(t/n 2 ) 

£ jjVSk„| S + ||l»ol + ^e^VV4 + („-V 2) . 

We now select so = An^ 1 ^ for some large constant A; this gives 

An-^ - E S0 Z- t <{\ + % + io^ M4 ) -~ 1/4 + 0(^ 1/2 ) • 

Choosing A large enough to swallow the constant C2 as well as the term 5 2 n 
(using the fact that 5 2 n is bounded), we obtain that 

E S0 Z- t > \au^ . 

Translating Zt back to \St\, we obtain 

E S0 |5 f | > rf ■ \atT x ' a > VAn~ l / A =: B , (5.4) 
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provided that A is sufficiently large (once again, using the fact that rf is 
bounded, this time from below). Since 

16 

Var^SV) < m/n 2 = , (5.5) 

the following concentration result on the stationary chain {St) will complete 
the proof: 

P^(|5 t | > An^ 4 )} < e{A) , and lim e{A) = . (5.6) 

A^oo 

Indeed, combining the above two statements, Chebyshev's inequality implies 
that 

||P So 0% G •) - vtHtv > tt([-B/2, B/2]) - P S0 (|^| < B/2) 

> i _ ^ _ e(x /Z) . (5.7) 

It remains to prove (5.6). Since we are proving a lower bound for the mixing- 
time, it suffices to consider a sub-sequence of the 8 n -s such that 5 n i/n con- 
verges to some constant (possibly 0). The following result establishes the 
limiting stationary distribution of the magnetization chain in this case. 

Theorem 5.1. Suppose that bim^oo 5 n yjn = a € R. The following holds: 

*{4-4) • < 5 - s > 

Proof. We need the following theorem: 

Theorem 5.2 ([8, Theorem 3.9]). Let p denote some probability measure, 
and let S n {p) = ^ X)j=i -^i(p)> where the {Xj{p) : j € [n]} /lave jomi 
distribution 

n 

Y[dp{xj) , 

i=l 

and Z n is a normalization constant. Suppose that {p n : n = 1,2,...} are 
measures satisfying 

exp{x 2 /2)dp n -> exp(x 2 /2)(i/9 . (5.9) 

Suppose further that p has the following properties: 

(1) Pure/ the function 

G p {s) :=j- log J e sx dp{x) 

has a unique global minimum. 

(2) Centered at m: let m denote the location of the above global mini- 
mum. 



exp 



{x 1 + . . . + x n ) 
2n 
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(3) Strength 5 and type k: the parameters k,5 > are such that 

G p (s) = G p {m) + 8 [S ~ 2 ™; + o((s - m) 2k ) , 

where the o{-)-term tends to as s — > m. 
If, for some real numbers a\, . . . , a2fc-i we have 



l-j/2fc 



+ o(n 



) , j = 1,2, ... ,2k — 1 , n ^ oo , 



i/ien i/ie following holds: 



and 



S n (Pn) ~ m 



n 



-l /2k 



S n {Pn) —> l{s^m} 
«1 1 



5 ' <5 



1 , 



exp 



. s 2 k 
\2k)\ 



2k-l 

£ 



O: 



ifk = \, 

if k > 2 . 



where 5 



1 > for k = 1 . 



Let p denote the two-point uniform measure on {—1,1}, and let p n denote 
the two-point uniform measure on {— n ,(3 n }. As |1 — (3 n \ = 8 n = 0(l/y/n), 
the convergence requirement (5.9) of the measures p n is clearly satisfied. We 
proceed to verify the properties of p: 



s 2 f s 2 s 4 

G P (s) = — — log / e sx dp(x) = — — logcosh(s) 



+ 0(s 8 ) . 



2 ° w 12 45 

This implies that G p has a unique global minimum at m = 0, type k = 2 
and strength 5 = 2. As 5 n y/n — ► a, we deduce that the G Pn -s satisfy 

s 2 

G Pn( s ) = -logcosh(/? n s) , 



G p ^0)=G^(0) = 0, 
G p 2 ](0) = l-(3 2 n = 5 n (2-5 n ) 



2a 



■n 



+ o{n~ 1 / 2 ) . 



This completes the verification of the conditions of the theorem, and we 
obtain that 

S n (Pn) 



exp 



„-l/4 

Recalling that, if Xi = ±1 is the i-th spin, 



s 4 s 2 

a — 

12 2 



, . . . , x n 



1 ( p ^ 



n 



(5.10) 
(5.11) 



l<i<j<n 
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clearly S n (p n ) has the same distribution as for any n. This completes 
the proof of the theorem. ■ 

Remark. One can verify that the above analysis of the mixing time in the 
critical window holds also for the censored dynamics (where the magne- 
tization is restricted to be non-negative, by flipping all spins whenever it 
becomes negative). Indeed, the upper bound immediately holds as the cen- 
sored dynamics is a function of the original Glauber dynamics. For the lower 
bound, notice that our argument tracked the absolute value of the magneti- 
zation chain, and hence can readily be applied to the censored case as-well. 
Altogether, the censored dynamics has a mixing time of order n 3 / 2 in the 
critical window 1 ± S for 5 = 0(l/y/n). 

5.3. Spectral gap analysis. The spectral gap bound in the critical tem- 
perature regime is obtained by combining the above analysis with results of 
[5] on birth- and-death chains. 

The lower bound on gap is a direct consequence of the fact that the mixing 
time has order ra 3 / 2 , and that the inequality i REL < i M ix(i) always holds. It 
remains to prove the matching bound t MDC Q) = 0(t REL ). Suppose that this 
is false, that is, t REh = o 

Let A be some large constant, and let sq = An^ 1 ^. Notice that the case 
5 2 n = 0(1) in Theorem 4.3 implies that Eito = 0(n 3//2 ). Furthermore, by 
Theorem 5.2, there exists a strictly positive function of A, s(A), such that 
lim^oo e(A) = and 

l -e(A) < tt(S > s ) < 2e(A) 

for sufficiently large n. Applying Lemma 3.7 with a = tt(S > sq) and 
(3 = h gives E SQ ro = o(n 3//2 ). As in Subsection 5.2, set t = n 3//2 /^4 4 for some 
large constant A. Combining Lemma 3.8 with Markov's inequality gives the 
following total variation bound for this birth-and-death chain: 

\\P S0 {S- t e •) " ttIItv < 4e(A) + o(l) . (5.12) 

However, the lower bound (5.7) obtained in Subsection 5.2 implies that: 

||P Sf) (S f e •) - 7t|| tv > 1 - 4e(VI/2) - U/A b . (5.13) 

Choosing a sufficiently large constant A, (5.12) and (5.13) together lead to 
a contradiction for large n. We conclude that gap = 0(n~ 3 / 2 ), completing 
the proof. 

Note that, as the condition gap • t Mlx (\) — > oo is necessary for cutoff 
in any family of ergodic reversible finite Markov chains (see, e.g., [5]), we 
immediately deduce that there is no cutoff in this regime. 
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Remark. It is worth noting that the order of the spectral gap at j3 c = 1 
follows from a simpler argument. Indeed, in that case, the upper bound on 
gap can alternatively be derived from its Dirichlet representation, similar 
to the argument that appeared in the proof of Proposition 3.9 (where we 
substitute the identity function, i.e., the sum of spins in the Dirichlet form). 
For this argument, one needs a lower bound for the variance of the stationary 
magnetization. Such a bound is known for /3 C = 1 (see [9]), rather than 
throughout the critical window. 

6. LOW TEMPERATURE REGIME 

In this section we prove Theorem 3, which establishes the order of the 
mixing time and the spectral gap in the super critical regime (where the 
mixing of the dynamics is exponentially slow and there is no cutoff). 

6.1. Exponential mixing. Recall that the normalized magnetization chain 
St is a birth-and-death chain on the space X = {—1, — 1+^, . . . , 1— ^, 1}, and 
for simplicity, assume throughout the proof that n is even (this is convenient 
since in this case we can refer to the state. Whenever n is odd, the same 
proof holds by letting ~ take the role of the state) . 
The following notation will be useful. We define 

X[a, b] := {x G X : a < x < b} , 

and similarly define X(a,b), etc. accordingly. For all x E X, let p x ,q x ,h x 
denote the transition probabilities of St to the right, to the left and to itself 
from the state x, that is: 

l-x 1 +tanh(/3(x + n -1 )) 



Px ■= Pm (x,x+ -) 
Qx ■■= Pm (x,x- I) 



2 2 

l+x 1 - tanh(/3(x - n' 1 )) 



2 2 
h x ■■= Pm (x, x) = 1 - p x - q x . 

By well known results on birth-and-death chains (see, e.g., [12]), the resis- 
tance r x and conductance c x of the edge (x, x + 2/n), and the conductance 
c' x of the self-loop of vertex x for x £ X[0, 1] are (the negative parts can be 
obtained immediately by symmetry) 

= ] [ — , C x = I I — , C x = ^ (C x _ 2 /n + Cx) , (6.1) 

ye*(0,x] Pv yeX(0,x] Qy Px + Qx 

and the commute-time between x and y, C X}V for x < y (the minimal time 
it takes the chain, starting from x, to hit y then return to x) satisfies 

ECs.j, = 2c s R(x <-> y) , (6.2) 
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where 

c s := ^2(c x + c x ) and R(x <-> y) := ^ . 

Our first goal is to estimate the expected commute time between and £. 
This is incorporated in the next lemma. 

Lemma 6.1. The expected commute time between and £ has order 

t exp :=^exp(|£log^||)dx, (6.3) 

where g{x) := (tanh(/3x) — x) / (1 — xtanh(/3x)). In particular, in the spe- 
cial case 5 — > we /iai>e ECb,£ = f exp ((| + o(l))5 2 n), where the o(l)-term 
tends to as n — > oo. 

Remark. If £ instead we simply choose a state in ^ which is the 
nearest possible to £. For a sufficiently large n, such a negligible adjustment 
would keep our calculations and arguments in tact. For the convenience of 
notation, let £ denote the mentioned state in this case as well. 

To prove Lemma 6.1, we need the following two lemmas, which establish 
the order of the total conductance and effective resistance respectively. 

Lemma 6.2. The total conductance satisfies 

<* = 6 ( A /?exp / ? log(i±44W 



8-*\2j ~»\l-g(x). 

Lemma 6.3. The effective resistance between and Q satisfies 

R(0 «-►£) = &(y/n/5) ■ 

Proof of Lemma 6.2. Notice that for any x £ X, the holding probability 
h x is uniformly bounded from below, and thus c' x can be uniformly bounded 
from above by [c x + c x _2/ n )- It therefore follows that cs = @(c~s) where 
cs := J2xex °xi an d it remains to determine %. We first locate the maximal 
edge conductance and determine its order, by means of classical analysis. 



logc - V log^- V logf 1 -^ l+tanh^ + n-im 

= E lo g (^ + 0(l/n))= £ logf^W) 

(6.4) 

Note that g{x) has a unique positive root at x = £, and satisfies g{x) > 
for x G (0, C) and (7(2;) < for x > (. Therefore, 

log c x < log c ? + O(x) < log c f + 0(1) , 
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thus we move to estimate cq. As logc^ is simply a Riemann sum (after an 
appropriate rescaling), we deduce that 



and therefore 



c f = 9 ( exp ( - 



n 



l + g(x) 
g(x) 



dx 



(6.5) 
(6.6) 



c x = 0(c c ) . 

Next, consider the ratio c x+2 / n /c x ; whenever x < £, g(x) > 0, hence we have 

Cx+2/n _ Px+2/n 1 + #(z) 



- 0(1 /n) > 1 + 2g(x) - 0(1/ n) . 



Whenever -^== < x < C — (using the Taylor expansions around and 



around C) we obtain that tanh(/3x) — x > \ y/S/n. Combining this with the 
fact that x tanh(/3x) is always non-negative, we obtain that for any such x, 
2<?(x) > y/S/n. Therefore, setting 




, 6 := C - 




(6.7) 



we get 



(6.8) 



> 1 + J- ~ 0(l/n) for any x £ . 
c x V fi 

Using the fact that 5 2 n — > oo, the sum of c^-s in the above range is at most 
the sum of £ 
position c^: 



the sum of a geometric series with a quotient 1 / '(1 + '$ / 'n) and an initial 



(6.9) 



We now treat x > £3; since g(() = and g(x) is decreasing for any x > (, 
then in particular whenever ( + y/5/n < x < 1 we have —1 = g(l) < g(x) < 
0, and therefore 

^±^ = ^±^<l + ff (x) + 0(l/n) . 

c x Qx+2/n 

Furthermore, for any ( + yfbjn < x < 1 (using Taylor expansion around Q 
we have tanh(/3x) — x < —yfsjn, and hence g(x) < —y/5/n. We deduce 
that 

c a;+2/n 



< 1 - \ - + 0(l/n) for any x e 1] , 
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and therefore 

Y cx<2i/|-cf. (6.10) 

XEX[&,1] 

Combining (6.9) and (6.10) together, and recalling (6.6), we obtain that 

c s = 2 (l x % &] I + 5v^+ I ^[6,6]|)cc = o(J^<^\ . 

Finally, consider x £ A? [£2 ,£3]; an argument similar to the ones above (i.e., 
perform Taylor expansion around £ and bound the ratio of c x+2 / n /c x ) shows 
that c x is of order in this region. This implies that for some constant b > 

c s > Y c - > fo l^[6,6]|c c >6W^c c , (6.11) 
xex[&,&] 

and altogether, plugging in (6.5), we get 

*- 8 (VMsjM£$)*))- (6l2) 



Proof of Lemma 6.3. Translating the conductances, as given in (6.8), to 
resistances, we get 

^t 2 ^ < 1 - \ - ~ 0(l/n) for any x G , 
r x \ n 

and hence 

^ < r^2y/n/5 < 2y/n/S , 

where in the last inequality we used the fact that r x < r x _ 2 / n (< r o = 1) for 
all x e X[0,(\, which holds since q x ^ p x for such x. Altogether, we have 
the following upper bound: 

r(o^o= Y r *+ Y r * + Y r * 

xeX[o,ti\ xexfafa] xeX[&,Q 

< |^[0,ei]| + 2w|+|Af[6,C]| <4v^- (6-13) 
For a lower bound, consider x G <Y[0,£i]. Clearly, for any x < we have 
= i-ftlnM/al) ^ 2fe ' and hence 

= 1 - g(x) + > 1-5x5 > exp(-6x<5) , 

r x 1 + g(x) 

yielding that 

r a > exp (-3<5n • £?) > e -3 . 
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Altogether, 

R(0 <-> C) > e~ 3 |^[0,Ci]| > , 
and combining this with (6.13) we deduce that R(0 *-* () = ®{y/n/5). ■ 

Proof of Lemma 6.1. Plugging in the estimates for &s and R(0 *-* () in 

(6.2), we get 

This completes the proof of the lemma 6.1. ■ 

Note that by symmetry, the expected hitting time from £ to — £ is exactly 
the expected commute time between and £. Hence, 

E c [r_ c ] = e(texp) • (6.15) 

In order to show that the above hitting time is the leading order term in 
the mixing-time at low temperatures, we need the following lemma, which 
addresses the order of the hitting time from 1 to £. 

Lemma 6.4. The normalized magnetization chain St in the low temperature 
regimes satisfies Eit^ = o(i ex p) • 

Proof. First consider the case where 5 is bounded below by some constant. 
Notice that, as p x < q x for all x > £, in this region St is a supermartingale. 
Therefore, Lemma 3.5 (or simply standard results on the simple random 
walk, which dominates our chain in this case) implies that Eir^ = 0(n 2 ). 
Combining this with the fact that t e xp > exp(cn) for some constant c in this 
case, we immediately obtain that Eir^ = o(i e xp)- 

Next, assume that 5 = o(l). Note that in this case, the Taylor expansion 
tanh(/3x) = j3x - l{f3x) 3 + 0((f)xf ) implies that 

C = ^35/ (3 3 - 0((/3C)) 5 = V35 + 0(5 3/2 ) . (6.16) 

Recalling that E[5t+i | St = s] < s + ^(tanh(/3s) — s) (as s > 0), Jensen's 
inequality (using the concavity of the Hyperbolic tangent) gives 

E[S t+1 - S t ] = E(E[S t+1 - S t | S t ]) < - (Etanh(/?5 t ) - E5 t ) 

n 

< - (tanh(/3E5 t ) - ES t ) . (6.17) 

n 

Further note that the function tanh(/?s) has the following Taylor expansion 
around C (for some £ between s and £): 

tanh(/3 S ) = ( + (3(1 - ( 2 )(s - C) + /3 2 (-l + C 2 )C(* - C) 2 

+ y (-1 + 4 C 2 " C 4 )(* - C) 3 + - C) 4 . (6.18) 
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Since tanh^ 4 ^(x) < 5 for any x > 0, (6.18) implies that for a sufficiently 
large n the term — g(s — () 3 absorbs the last term in the expansion (6.18). 
Together with (6.16), we obtain that 

tanh(/fc) < C + P(l ~ ( 2 )(s - C) + (3\-l + C 2 )^(s - C) 2 • 
Therefore, (6.17) follows: 

HSt+i - S t ] < ~(ES t - C) 2 . (6.19) 
In 

Set 

bi = 2~ % , %2 = min{i : bi < Vd} and u% = min{t : ESt — ( < bi] , 

noting that this gives b{/2 < E5( — ( < bi for any t G [ui,Ui + \\. It follows 
that 

bJ2 An 
Ui+l ~Ui < — = —=- , 

2n V 2 I % 

and hence 

^ - Uj < ^ = 0(n/tf) , 

i=l i:bf>8 1 

where we used the fact that the series {b^ 1 } is geometric with ratio 2. We 
claim that this implies the required bound on EiT£. To see this, recall (6.19), 
according to which Wt ■= n(St— C)l{ Tc >t} is a supermartingale with bounded 
increments, whose variance is uniformly bounded from below on the event 
T£ > t (as the holding probabilities of (St) are uniformly bounded from 
above, see (3.4)). Moreover, the above argument gives EW t < nVS for some 
t = 0(n/5). Thus, applying Lemma 3.5 and taking expectation, we deduce 
that EiT^ = 0(n/ 8+5n 2 ) = 0(5n 2 ), which in turns gives Eir^ = o(t exp )- ■ 

Remark. With additional effort, we can establish that Eot+£ = o(t ex p) (for 
more details, see the companion paper [6]), where t±^ = min{i : | jS* | > C}- 
By combining this with the of St symmetry and applying the geometric trial 
method, we can obtain the expected commute time between and £: 

E C T = (i + o(l))EC , c = e(t cxp ) , 

and therefore conclude that EiTq = &(t exp ). 
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6.1.1. Upper bound for mixing. Combining Lemma 6.4 and (6.15), we con- 
clude that Eit_£ = 0(t e xp) and hence EiTo = 0(t exp ). Together with 
Lemma 3.1, this implies that the magnetization chain will coalescence in 
O(texp) steps with probability arbitrarily close to 1. At this point, Lemma 
3.3 immediately gives that the Glauber dynamics achieves full mixing within 
0{n log n) additional steps. The following simple lemma thus completes the 
proof of the upper bound for the mixing time. 

Lemma 6.5. Let t exp be as defined in Lemma 6.1. Then nlogn = o{t exp ). 

Proof. In case 5 > c > for some constant c, we have t exp > nexp(c'n) for 
some constant d > and hence nlogn = o(t exp ). It remains to treat the 
case 5 = o(l). 

Suppose first that 5 = o(l) and 5 > cn -1 / 3 for some constant c > 0. In this 
case, we have t 0X p = f exp ((| + o(l))<5 2 n) and thus nexp(|n 1 / 3 ) = 0(t exp ), 
giving nlogn = o(t exp ). Finally, if 5 = o(n -1 / 3 ), we can simply conclude 
that n 4//3 = 0(i ex p) and hence nlogn = o(t exp ). ■ 

6.1.2. Lower bound for mixing. The lower bound will follow from showing 
that the probability of hitting — £ within et exp steps is small, for some small 
e > to be chosen later. To this end, we need the following simple lemma: 

Lemma 6.6. Let X denote a Markov chain over some finite state space Q, 
y £ fi denote a target state, and T be an integer. Further let x £ O denote 
the state with the smallest probability of hitting y after at most T steps, i.e., 
x minimizes P 1! (r !/ < T). The following holds: 



Proof. Set p = P x (t v < T). By definition, P z (r y < T) > p for all z G ft, 
hence the hitting time from x to y is stochastically dominated by a geometric 
random variable with success probability p, multiplied by T. That is, we 
have E x T y < T/p, completing the proof. ■ 

The final fact we would require is that the stationary probability of 
X[—l, — C] is strictly positive. This is stated by the following lemma. 

Lemma 6.7. There exists some absolute constant < C n < 1 such that 



CV<7r(*[C,l]) (=7r(*[-l,-CD) • 



Proof. Repeating the derivation of (6.11), we can easily get 




) 



dx 



) 
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Combining the above bound with (6.12), we conclude that there exists some 
CV > 0, such that n(X[(, 1]) > C v . ■ 

Plugging in the target state — £ into Lemma 6.6, and recalling that the 
monotone-coupling implies that, for any T, the initial state so = 1 has the 
smallest probability (among all initial states) of hitting —( within T steps, 
we deduce that, for a sufficiently small e > 0, 

Pi(t_ c < rf cxp ) < \c n . 



2 



This implies that 
which in turn gives 



^exp 



O mix ^ — Cjt j J , 



^exp — O (^mix^)) 



6.2. Spectral gap analysis. The lower bound is straightforward (as the 
relaxation time is always at most the mixing time) and we turn to prove 
the upper bound. Note that, by Lemma 6.7, we have tt(X[£, 1]) > C n > 0. 
Suppose first that gap • i M ix(z) — * °°- I n this case, one can apply Lemma 
3.7 onto the birth- and-death chain (St), with a choice of a = ir(X[£, 1]) and 
(3=1- tt(X[(, 1]) (recall that W(|) = 0(Eit_ c )). It follows that 

E c r_ c = o (Eir-f ) • 

However, as both quantities above should have the same order as i M ix( \), this 
leads to a contradiction. We therefore have gap-£ MIX (^) = 0(1), completing 
the proof of the upper bound. 

References 

[1] D. Aldous and J. A. Fill, Reversible Markov Chains and Random Walks on Graphs. 

In preparation, http: //www. stat .berkeley . edu/~aldous/RWG/book .html. 
[2] M. Aizenman and R. Holley, Rapid convergence to equilibrium of stochastic Ising 

models in the Dobrushin Shlosman regime, (Minneapolis, Minn., 1984), IMA Vol. 

Math. Appl., vol. 8, Springer, New York, 1987, pp. 1-11. 
[3] R. Bubley and M. Dyer, Path coupling: A technique for proving rapid mixing in 

Markov chains, Proceedings of the 38th Annual Symposium on Foundations of Com- 
puter Science (FOCS), 1997, pp. 223-231. 
[4] M.-F. Chen, Trilogy of couplings and general formulas for lower bound of spectral 

gap, Probability towards 2000 (New York, 1995), Lecture Notes in Statist., vol. 128, 

Springer, New York, 1998, pp. 123-136 1157,230,231. 
[5] J. Ding, E. Lubetzky, and Y. Peres, Total-variation cutoff in birth- and- death chains, 

preprint. 

[6] J. Ding, E. Lubetzky, and Y. Peres, Censored Glauber dynamics for the mean-field 
Ising Model, preprint. 

[7] R. S. Ellis, Entropy, large deviations, and statistical mechanics, Grundlehren der 
Mathematischen Wissenschaften, vol. 271, Springer- Verlag, New York, 1985. 



MIXING TIME EVOLUTION OF GLAUBER DYNAMICS 



43 



[8] R. S. Ellis and C. M. Newman, Limit theorems for sums of dependent random variables 
occurring in statistical mechanics, Z. Wahrsch. Verw. Gebiete 44 (1978), no. 2, 117- 
139. 

[9] R. S. Ellis, C. M. Newman, and J. S. Rosen, Limit theorems for sums of dependent 
random variables occurring in statistical mechanics. II. Conditioning, multiple phases, 
and metastability, Z. Wahrsch. Verw. Gebiete 51 (1980), no. 2, 153-169. 

[10] R. B. Griffiths, C.-Y. Weng, and J. S. Langer, Relaxation Times for Metastable States 
in the Mean-Field Model of a Ferromagnet, Phys. Rev. 149 (1966), 301 -305. 

[11] D. A. Levin, M. Luczak, and Y. Peres, Glauber dynamics for the Mean-field Ising 
Model: cut-off, critical power law, and metastability, preprint. 

[12] D. Levin, Y. Peres, and E. Wilmer, Markov Chains and Mixing Times, 2007. In 
preparation, available at http://www.uoregon.edu/~dlevin/MARKOV/. 

[13] §. Nacu, Glauber Dynamics on the Cycle is Monotone, Probability Theory and Re- 
lated Fields 127 (2003), 177-185. 

Jian Ding 

Department of Statistics, UC Berkeley, Berkeley, CA 94720, USA. 
E-mail address: jding@stat.berkeley.edu 

Eyal Lubetzky 

Microsoft Research, One Microsoft Way, Redmond, WA 98052-6399, USA. 
E-mail address: eyal@microsoft.com 

Yuval Peres 

Microsoft Research, One Microsoft Way, Redmond, WA 98052-6399, USA. 
E-mail address: peres@microsoft.com 



